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ABSTRACT 

Self similar 3D distributions of point-particles, with a given quasifractal di- 
mension were generated on a Menger sponge model and then compared with 
2dfGRS and Virgo project data Q. Using the principle of local knowledge, it is 
argued that in a finite volume of space only the two-point minus estimator is 
acceptable in the correlation analysis of self similar spatial distributions. In this 
sense, we have simplified the Pietronero-Labini correlative analysis by defining a 
/T-minus estimator, which when applied to 2dfGRS data revealed the quasifrac- 
tal dimension D ^ 2 , as expected. In our approach the fC-minus estimator is 
used only locally. Dimensions between D = 1 and D = 1.7, as suggested by 
the standard ^(r) analysis, were found to be fallacy of the method. In order to 
visualize spatial quasifractal objects, we created a small software program called 
RoPo ("Rotate Points"). This program illustrates and manifests local correlative 
analysis in which the visual inspection emerged as a first step and a key part of 
the method. Finally, we discuss importance and perspective of the visual inspec- 
tion on available real and simulated distributions. It is also argued that results 
of contemporary cosmological simulations do not faithfully represent real data, 
as they show a formation of ever increasing coUapsars. We consent that 2dfGRS 
data are reminiscent of some kind of underlying turbulence like effects in action. 

Subject headings: 3D random fractals, Menger sponge, self-similarity, quasifractal 
dimension, principle of local knowledge, correlation analysis, minus estimator, 
2dfGRS, Virgo project simulations, visual inspection. 



www . mso . anu .edu.au/2dFGRS/, www . mpa-gar ching. mpg . de / Virgo / 
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Introduction 



A statistical app roach to the analysis of galaxy and cluster correlations (like in 



Gabrielli et al. 



(120051 )) has led to the surprising result that they may be fractal up to the 
limits of the available catalogs. This so called fractal conjecture is based on the analysis 
of two-point correlation function, presently the main statistical tool in the observational 
cosmology. From this point of view and together with the emersion of modern redshift 
surveys, the galaxy structures appear highly irregular and self-similar with fractal dimension 



D ~ 2 up to the deepest scales probed so far (lOOO/i ^Mpc 
indicated from the interpretation of the number counts as in 



and maybe even more as 



Labini et al 



fll998h . 



Clearly and obviously, space is highly structured. However, a meaningful explanation 
for this large scale structure (LSS) is still lacking, and this remains the key question in 
contemporary cosmology. Its resolution seems crucial for possible extension of general 
relativity to the quantum domain and so on. 

In the following, respecting the fractal conjecture, we comprehend and repeat the 
correlation analysis of available data. We use the i^-minus estimator as a tool in spatial 
statistics when dealing with selfsimilar structures. 

As usual, the three-dimensional distribution of galaxies and generated quasifractal 
objects are represented by a set of points in empty space. Their density is given by the sum 
of delta functions n{r) = ^ 5(r — ri) so reducing all properties of these objects to simple 

i 

material (mathematical) points, defined only by their positions in space. Of course, it is 
always possible to add multiple attributes (mass, "colour", orientation, ...) to these points 



a 



^In a perspective even thousands of simple or extended attributes can sit in the back- 
ground, organized in some database model. Variety of metrics can be conceived, etc. Con- 
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Our starting relation is N{r) ~ r^, known as a "mass-length" relation (IBaveye et al. 



( Il998l )). We imagine a sphere of radius r around every point of a quasifractal object and 



count all the points inside the sphere to get 



iV(r) = 5^ ^^(r - = Ntrendir) + Nfiucir) 



The "mass-length" relation refers to the behaviour of Ntrend- In further discussion ~ will 
be used in a sense of trend, or mean. 

A quasifractal object is made up of a large but finite number of points. It is usually 
only approximately self-similar at some scale range. The scale is necessarily restricted 



from above by t he size of the 



more details see 



Baveye et al. 



s amp 



e, and from below by the size of the galaxy itself (for 



(119981 ) ). Such a scale bounded object will be referred to as 



quasifractal to distinct it from an exact fractal which is composed of an infinite number of 
points and is exactly self similar at all scales, in other words, it is infinitely large in extend 
and infinitely structured [f]. The "mass-length" relation of a quasifractal object has a power 
law trend only in some range of scales. 

An interesting question arises now, how these random three-dimensional quasifractal 
objects really look like. To answer this question, we created a small program for generating 



various such objects using the Menger sp onge model (MSM) 



determined by the formula Nr^ = 1, as in 



Baveye et al 



The fractal dimension is 



(Il998l ) so that we need only to give 



nections also may have attributes, we will call them "predicates" or "relations", but this is 
another subject not to be discussed now. 
^reader may try GNU XaoS software 

^program is named "menger. for" and is available upon request from martims@irb.hr ; it 
does 5 or 6 selfsimilar recursive iterations depending either on the choice of N and 1/r, or 
on the available number of sample points 



- 5 - 



two numbers, and 1/r, respectively. 

By using an uniform random number generator for position coordinates, we also 
generated an "uniform gas" of particles. In our discussion, it is regarded as a quasifractal 
object too, with fractal dimension D = 3. Although an uniform gas is a trivial quasifractal 
object, it is important in the correlation analysis for it has no imposed or presumed inner 
structure. As such, it is suitable for basic testing of various two point estimators. 



Our correlation analysis is bas ed on the siraplifica tion of the Pietronero-Labini r(r 



and r*(r) functions as presented in 



Labini et al. 



(119971): 



r+Ar 



rir) 



N ^ Anr'^Ar 

i=l 



n{fi + r')dr' ^ 



r*(r) 



We make use only of the integral kernel of the r(r) function parameterized by Ar and 
define our K function as: 



K{r, Ar) = ^ Kf^r, Ar) = ^ / n{fi + r')dr' ~ r^"^ 

and r < |r^| < r + Ar. This is our main formula. It can be easily generalized to 
the multidimensional case. i^'(r, Ar) and Kff^{r,Ar) are ~ (trend) in the case of 

quasifractal. 

Simply as it is, Ki{r, Ar) counts the number of points falling into the spherical shell 
(r, r + Ar) around some starting point r, (the origin). We sum over all such starting points 
that can be found in the volume of the sample. In fact K{r, Ar) counts the number of pairs. 
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or the number of connections between the points. Any correlation analysis basically count 
and classify connections between points. Obviously, K is the simplest way of doing that. 

Pair connections and their possible attributes are usually classified by length only, 
though we may conceive many other classification criteria. 

No additional statistical or probabilistic interpretation of K{r^ Ar) is needed, in 
particular there is no need for a normalization, sample defined multipliers, or for the 
problematic factors and 1/Ar as in r(r) and r*(r) functions, respectively. We are 
interested only in the trend of K, expected to be ~ r^~^ and in the determination of Z) - 
we quit any additional statistical interpretation of correlation analysis apart from classifying 
connections. This simplification is found to be very convenient in the correlation analysis 
of highly irregular selfsimilar spatial distributions. 

In most of the present literature it is debated how to define average densities of irregular 
distributions. Here, we make a philosophical detachment and debate about connections 
(pair links) only, avoiding density as much as possible. It is a well known fact, that an 
"average density" cannot be defined for a selfsimilar (fractal) distribution as it depends on 
the scale and the sample volume. In the case of a mathematical fractal, the average density 
is zero ("dominated by space") and the only exception is an uniform gas (D = 3). Therefore 
and furthermore, fluctuations are also a problem - fluctuations around what? There is 
neither a reference object nor a reference s cale - all scales are relevant. We generally avoid 



(119971 )). and anything derived from 



using terms like homogeneity and isotropy (ILabini et al. 
them. 

Estimator is a particular implementation of a correlation function - let's say numerical 
procedure. Examining formula (1), we recognize immediately the necessity that all spherical 
shells remain entirely in the volume of the sample. This means that some points in the 
sample near the edge cannot be used as starting points for a given shell (r, Ar) - there 
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are shells not fully in the sample. From this observation, we construct the so called minus 
estimator K{r,Ar). 

For a reasonably small Ar, the trend of the K{r, Ar) function behaves like a derivative 
of the original mass-length relation N{r) ~ r^. Of course, we didn't specify this case in the 
theory, but in practice this is the case when constructing histograms. 

We notice that all information r(r), r*(r) and ^(r) can offer are already fully 
contained in K{r,Ar). For example, ^{r) = 1, with r(r) as a key part, and 



-1 as an unimportant part ("noise"), has a trend ~ r 



D-3 



small compared to th e first part 



extended discussion in 



Baryshev et al. 



only in the region wh ere 



Baryshev fc Teerikorpil (120051 ) 



IS 



(119981 ): iBaryshevi (119991 )). Among many other flaws, ^(r) has also a problem 
with the definition of an average density Uaverage which cannot be properly defined in 
the case of quasifractal distributions. Nevertheless, ^(r) can be applied to continuous 
distributions with small (or very small) fluctuations of average density, such as dense fluids 
(or perhaps to cigarette smoke like distributions), but not to self-similar distributions. We 
wonder why majority of authors, in spite of problems, prefer to use ^(r). 



2. K{r, Ar) as a minus estimator and construction of histograms 

A spatial pattern of galaxies may be virtually infinite in extent, but we can observe 
only part of it through a restricted window ( "sample" ) which give rise to a sampling bias 
and censoring (missing data). When a spatial point distribution is observed through a 
bounded window, edge effects hamper both the information about the properties of the 
pair distance distribution and on the function K{r,Ar). In fact, the distance from a given 
starting point to the next point of the sample is censored by its distance to the boundary 
of the sample. Edge effects become obviously more severe as the distance r increases. 



- 8 - 



In a sample with points there are A^(A^ — l)/2 pair distances. From this set of pairs, 
we can use only a subset defined by pairs of points with at least one point which is more 
than d ( "maximum range" ) away from the nearest sample boundary. 

In spatial statistics, one traditionally uses the edge- corrected estimators which are 
weighted empirical distributions of the observed pair distances. The simplest approach 
is the "border method" , where attention is restricted to those starting points lying more 
than d units away from the boundary. All pair distances up to d are observed without 
censoring. This approach is often justified by appealing to the "local knowledge principle" 
of mathematical morphology stating that we can know only what is in the sample, and 
draw conclusions only from that knowledge. However, the border method usually discards 
much of data points, especially in three dimensions that can be unacceptably wasteful when 
estimating K. 



If we take into account the 
Baddeley et al. 



ocal knowledge principle, thoroughly discussed by 



Baddeley &: Gill! (jl997l ). we can avoid making any homogeneity or isotropy 
presumption. The use of i^'(r, Ar) as a minus estimator ("the border method estimator") 
seems unavoidable on quasifractal samples. Also, this estimator is very simple to implement. 
An observer must first make a shift d from the edge of the sample in order to choose the 
first point of the pair. The second point can be anywhere in the sample, but no more than 
d units away form the first point. 

The sub-sample formed by first points is, of course, smaller than the whole sample, 
and there is a compromise between the number of initially chosen points and the range d 
that can be varied. 

Edge effects are easily recognized on histograms, by decline of the graph, as soon as 
we approach the boundary of the sample (see Fig. 4.). There may be numerous histogram 
construction schemes, for example step of the histogram may be Ar or a fraction of it 
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for a smoother trend and so on. In calculations, we truncated decimals to get the index 
immediately, as in an appended procedure (see c source code in appendix) which is almost 
trivial, especially compared to other estimators. 

Generally, we avoid any edge corrections in treating apparently (or conjectured) 
selfsimilar or fundamentally irregular distributions. We simply don't know and cannot know 
what is outside of the sample volume where there may exist strong clustering of galaxies 
as well as big voids. In fact, only information that are contained in the sample volume 
should be used. Any estimator using weighting schemes or other kind of edge corrections 
necessarily puts in the calculation additional information that do not exist, and have no 
justification particularly in the case of random fractals or extreme clustering, making the 
calculating results more or less arbitrary. 

It will be later argumented that even if we use the minus estimators, as presently the 
best tools known to us, we cannot completely rely on the results of correlation analysis, 
which means that we cannot and should not draw any firm conclusions from these results, 
perhaps only to answer the question whether or not the observed distribution shows 
self-similarity in some scale range. This answer is given by reading histograms only, with 
no further conclusions. 

An example of a good application could be the mutual comparison of (slightly) different 
distributions. The correlation analysis yields a correlation histogram, which is a highly 
reduced and abstracted, imperfect and incomplete information about the distribution. Full 
information is, of course, given by specifying all points positions and their properties. In 
that sense, two different distributions may have the same histogram, but they are not 
even close in appearance. It is clear now that visual inspection must be also a key and 
unavoidable part of the correlation analysis. For this reason, we should first visually 
examine a given distribution, before any calculation, and see if it is appropriate for applying 
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K{r, Ar) which means that the chosen sample has no hollow in the middle where we search 
for first points, etc.. A histogram is fair (trend can be recognized) even for a small number 
of points. We call this "7V(7V - 1)" effect. 

A special and rare, unlike feature of a spatial point distribution is its self-similarity 
with a nontrivial quasifractal dimension D < 3 . If this feature exists, then there must be 
a mechanism or reason responsible for generating and supporting this fractal clustering 
(fractalisation with stable self-similarity) and not some other conceivable kind of clustering. 
Cases are just few: gas, collapse, fractal, then their combinations, and then something 
neither of these. 

For a self-similar distribution, we expect a histogram trend to follow the power-law 
behaviour. The standard straight line adjustment using the least square method in a log-log 
scaled diagram, gives according to the K{r, Ar) ~ j-^-^ formula, a power-law exponent, 
interpreted as D — 1. The adjustment to other functional forms may be also conceived, 
mainly to parabola, to distinguish coUapsars from a foam. 

The traditional correlation analysis uses all sample points equally and simultaneously, 
and in this sense, it is "global". Our approach is basically "local", it implies the use 
of sub-samples, even edges needn't to be exactly defined, although we should take care 
of the approximate range d. Any general conclusion about the whole sample should be 
avoided. Instead, we examine a correlative histogram around each point fi determined 
by Kff^{r,Ar) and casually parameterized by Ar. We know from experience that this 
histogram is changing quickly from one point to another - so that giving too much weight 
to the correlative analysis is not necessary nor recommended. As an illustration our local 
correlative analysis is implemented in the RoPo program together with the possibility of 
local sampling. 

In most cases we can choose one maximal (cubical) sub-sample, and then discuss 
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correlative analysis of the sample, not the analysis around each point. The calculation is 
demanding because of the great number of connecting pairs, so in order to limit calculations 
to a reasonable time span of few seconds, we used the Monte Carlo sampling procedure. It 
has not been proved yet, but it appears that a preliminary trend can be obtained already 
within the first (mill) seconds. 

The advantages of K{r, Ar) minus estimator are the following: the histogram is always 
positive, it has a rising trend like parabolas (~ r'^) , much easier perceived and interpreted 
than hyperbolas, which require log-log scaled diagrams. Further advantages are: fast and 
simple calculations compared to other estimators, and small hardware requirements. There 
is no need for additional distributions, normalizations, weighting schemes, etc. We believe 
it is the simplest method for testing self-similarity in spatial distributions of points. 

A prona ii ient and oft en mention case, relevant in cosmology, following for example 



Labini et al 



(11997 



expected power-law behaviour ~ r'' is a straight line z/ = 1 , 
parabola u = 2 (meaning D = 3), or something in between!^ 



19981 ). and others, is a quasifractal dimension D = 2 . On such 
distributions, our estimator gives a straight line, not only in log-log scaled diagram but 
also in the histogram itself. The interpretation of a histogram is very simple; we ask if the 

meaning D = 2), or a square 
Other cases are considered as 
deviations. Almost always a simple fast and rough visual estimate is sufficient to determine 
approximate quasifractal dimension D. It will be explained later, why more than that 
shouldn't be even asked. 

A little bit of philosophy; to us scaling exists only after the correlation analysis and 
histogram inspection have been completed on a chosen sample. This depends obviously on 
the sample size and its properties. Then, if we take another similar sample from the same 



< 1.8 means a rather extreme clustering, we don't expect to find nothing such. 



- 12 - 



distribution, and find no scaling, the question is whether the scahng exists or not? The 
only reasonable answer would be that locally on a particular sample, using a particular 
estimator, we see the scaling property. 

On a gas- like distribution K{r,Ar) is almost exact (square parabola), and this is the 
only case. Therefore, we can use K{r, Ar) to determine whether or not and how much a 
certain spatial distribution is gas-like (homogenous), and this is probably the best use of K. 

3. Compact samples and local analysis 

Prom a given large sample we pick up several sub-samples according to the following 
criteria. In order to avoid homogeneity and isotropy presumptions (anything with 
orientation should be carefully avoided) , we require that a chosen sub-sample looks compact 
in its visual appearance. A criteria of minimal surface area to volume ratio can be 
established. An elongated sub-sample, a filament for example, has obviously a well defined 
orientation in space and it can fall along and/or across the sample creating the problem 
for the estimator. It is rejected in our analysis because it calls for an isotropy presumption 
(equivalence of orientations). A spherical sample is the best example of a compact sample 
we are talking about. Prom the point of the calculation cost and speed, wc prefer use of 
cubically shaped samples. They do have a defined orientation, but in practice this is found 
to be unimportant. 

Elongated distributions such as 2dfGRS, the best available data, and SDSS, were 
analysed by first choosing a number of cubical samples, which may or may not intersect 
each other, as it makes no difference in the final results. Then we found a trend for each 
of the cubic samples and took an average value of them, eventually rejecting some of the 



- 13 - 

results that were too discursive for some reason 1^. Initial samples were not of the same size 
(a greater redshift 2; ~ r, bigger sample), but we rescaled them to the unit cube. As we 
presumed self-similarity of 2dfGRS data, this rescaling made no difference, the structure 
should be the same on all scales. In this way, the points (data) deficiency at larger z has 
also been avoided. We didn't try to correct any flaws of available data, nor did we select 
any particular distinguishing sub-set from them, by any criteria. It appears that correcting 
flaws of data will not bring any significant change in the final results. 



4. Special cases and application 

Let us now make few comments on our analysis of generated self-similar distributions. 
In the case of an uniform gas, the error of the calculated quasifractal dimension was in 
the second decimal digit, or higher, while in the case of a quasifractal object with D = 2, 
the error was in the first decimal digit. In order to obtain the result with an error in 
the second decimal digit, we had to generate and analyse about ten and more equivalent 
distributions generated with different random number seeds, and then take the average 
to discern imposed quasifractal dimension. Smaller quasifractal dimensions always mean 
stronger (more extreme) clustering for a given number of points in the sample and effects of 
the finite sample size and edge effects become more and more important. A lack of available 
starting points is also a problem. As we move away from D = 3, any kind of estimator 
becomes less and less reliable. For D = 3 all estimators may be used but for D < 2.6 
(estimated) we suggest only minus estimators. If there is any sense to talk about the error, 
it will be -D = 2.0 ± 0.3 in comparison with gas where we have D = 2.988 ± 0.064, as tried 
on examples. 

^software Korela has been made for processing a large number of samples in a sequence 
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From the above discussion, it has become clear that in a correlative histogram we 
cannot rely on one sample only and some kind of averaging over several equivalent samples 
is needed. For example, we can divide a big sample to several smaller subsamples. There 
are many possibilities. A series of evenly or randomly oriented samples can be chosen. Such 
calculation has been conducted on 2dfGRS data, to extract 259 samples on the first branch 
(some rejected) finding average D = 2.016, and 195 samples on the second branch finding 
average D = 1.976. These findings are consistent with D ^ 2.0 found by other authors (for 
example Pietronero and Labini). 

We are not inclined to treat this result as a final truth, but only as an illustration. 
Furthermore, the range of correlation analysis is in this case quite small, we need much 
bigger compact samples. From the point of view of our method, we can only draw 
conclusions up to this small range. 

The large-scale structure (LSS) of the universe has been heavily debated in the recent 
history without a firm conclusion. We believe that LSS, as well as interstellar g 

n 

represented by random fractals. From 2dfGRS data, LSS is obviously not homogeneous Q 
The question still remains whether LSS is self-similar or not. Some th eoretical aspects of 



possible self-similar LSS are discussed in 



Baryshev fc Teerikorpil (120051 ) 



''Upon request RoPo freeware will be given; Win32 for now, gcc port also exists. Software 
has been made for (serial) visualization of points distributions (simulation snapshots, or just 
a general point distributions, attributes also may be added - colour for example). 2df.b 
data will be provided in a RoPo .h format. Samples as big as 10 million points have been 
successfully visualized on a standard desktop computer. Main design objective was, of course, 
speed. 
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Conclusion 

As big compact samples cannot be extracted from 2dfGRS, only the short range 
i^'-minus estimator correlation analysis is possible. A visual inspection remains the only 
alternative on a greater scale, and we think there is self-similarity on the entire available 
scale range and there is no visible crossing to homogeneity within 2dfGRS {z = 0.26). We 
didn't treat less important effects like the F ingers-of-God, defi ciency of points on greater 



z, errors in z measurements, and so on, see I Jones et al. 



20051 ). A correlation estimator 



is a rough tool, as the reader may have noticed from our previous discussion, so we don't 
consider this to be an important issue. It is the same with the visual inspection method - 
flaws of data may be ignored for the first try. 

Visual comparison with generated quasifractal objects and simulations data has been 
also conducted. From it we decided by consensus that 2dfGRS data resembles fractal 
dimension between D = 2.1 and 2.2. What we see is similar to a foam (or maybe "drunk 
spider" net), and not to various Virgo project simulation distributions which are dominated 
by collapsars of different size. If those simulations do not provide something like 2dfGRS 
data, it means there is no explanation for this structure in the contemporary computational 
cosmology. There is nothing new in the theory either, although many alternatives and 
combinations of them have been suggested, MOND for example. From our perspective, 
some kind of turbulence should be incorporated in the simulation procedure (see picture 
of Crab Nebula's peripheral regions as an example). Of course, GTR should be a natural 
starting point. 

Visual inspection is proposed here as an important part of the correlation analysis, but 
also as a sole method. This involves psychology and the knowledge on perception of density 
contrast seen on the screen, in our case bright and dark areas or groups of pixels on the 
black background. Perception is of course nonlinear. In addition, pixels may overlay, etc.. 
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so we don't have definite and precise answer to those issues I^. But for the mere comparison 
of 2dfGRS and simulated data RoPo is sufficient, and we have firm conclusion that these 
distributions are different. 

The entire discussion about possible scale of transition to homogeneity should be 
postponed until at least one big compact sample became available, and for that purpose 
2dfGRS equipment could be reused. We think that present data do not show any transition 
to homogeneity. The fractal universe seems to be an alternative. 

Why did we create the RoPo program? A simple reason was that the human brain 
and eyes perceive objects in motion much easier and deeper than just plain static pictures. 
We tried to emulate advantage of stereoscopic seeing (standing eyes but moving objects 
- instead of standing objects but moving eyes as in nature), and it made a definite 
breakthrough in the discussion of 2dfGRS data morphology against the Virgo project 
(gravodinamical) simulation data morphology. In these simulations we clearly see two 
distinct classes of structures: coUapsars and filaments. What is even more interesting, the 
filaments are almost all straight line alike, what we recognize as intersections of planes 
caused by gravity amplified inflicted plain wave perturbations. Those distributions are 
utterly dominated by collapsars, contrary to 2dfGRS where net-like structures dominate. 

To resolve these problems more decidedly, we obviously need higher order correlation 
analysis. However at the moment there is no convincing way of defining three-points and 
other multi-points correlation functions in K{r, Ar) paradigm. The question is whether 
this attempt would be worth doing it, because the two-point correlative analysis is already 
too approximate and can not produce a firm conclusion. It gives only an approximate and 

^interesting subject, for example in automated medical imaging analysis, or maybe com- 
puter game industry 
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unreliable abstraction of properties of an underlying point process. 
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tO=t[i] [0] ;if (t0<0.25| | tO>0 . 75)goto jump; 

tl=t[i] [1] ;if (tl<0.25| | tl>0 . 75)goto jump; 

t2=t[i] [2] ;if (t2<0.25| 1 12>0 . 75)goto jump; 

/* second point */ j=(long) (n*genrand_real2() ) ; 

kO=tO-t [j] [0] ;kl=tl-t [j] [1] ;k2=t2-t[j] [2] ; 

k= (long) (1000. *sqrt (k0*k0+kl*kl+k2*k2) ) ; // index 

if ((k>250) I I (k<5))goto jump; ul [k] =ul [k] +1 ;nl=nl+l ; jump:;} 
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Fig. 1. — 3D Quasifractal object with fractal dimension D = 1.28 containing 46656 points, 
as generated by a Menger sponge model. Distribution of points is dominated by large voids 
- clustering is extremely strong. A lot of space had to be searched to find new points; see 
"mass-length" relation. In the older literature one of the usual errors was comparing this 
with 2dfGRS data. Obviously they are not equal. 



Fig. 2. — Quasifractal objects with fractal dimensions (a) D = 2.00, (b) D = 2.52 and (c) 
D — 3.0, respectively. Each picture has 1048576 points. We emphasize this number so that 
a reader can get impression of how self-similar clustering look like. All distributions are 
three-dimensional (note the shading) and they fit into the unit cube (menger generate by 
recursive and (quasi)random "slicing" from the unit cube). Distributions are rotated for a 
better view. 
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Fig. 3. — Correlation histogram of an uniform gas in a unit cube, taking into account the 
edge shift (lower straight line, D — 3.00, histogram accumulated less points), and without 
edge shift (upper curve, straight line adjustment gives D — 2.85 - obviously wrong). Notice 
that the correlation range is only 0.25. A httle difference in view, but great in numbers. 
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Fig. 4. — Behavior of the i^-cstimator on a Virgo project simulation sample - the trend 
resembles a parabola, and the distribution histogram is gas-alike, at least locally around an 
initial point - this distribution looks almost self-similar, but with D 2.63, not D 2 as 
expected; we also sec the edge influence on the right hand side. Using RoPo, we can move 
fast through a given distribution from one point to another, and on each position instantly 
find a histogram. 
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Fig. 5. — Clustering in 2dfGRS data; looks like points are running away from voids and this 
is a discrepancy with simulations data where points tend to collapse - cores can be clearly 
distinguished, see next figure. 
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Fig. 6. — Simulations - by correlation analysis more like a gas than 2dfGRS. If we replace 
clusters with points - they are even more like a gas. 
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Fig. 7. — SDSS data flaw - red shifts are rounded at a third decimal point and as we can 
take only small samples we observe slices - estimator clearly recognizes this artificial effect 
giving D pa 1.8 instead of expected D k,2. 




Fig. 8. — Menger sponge - picture taken from www.sewanee.edu/physics/ PHYSICS123/ 
physicsl23.html 



